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(Dated: September 25, 2007) 

We present a measurement of the muon charge asymmetry from W boson decays using 0.3 fb _1 
of data collected at y/s = 1.96 GeV between 2002 and 2004 with the DO detector at the Fermilab 
Tevatron pp Collider. We compare our findings with expectations from next-to-leading-order cal- 
culations performed using the CTEQ6.1M and MRST04 NLO parton distribution functions. Our 
findings can be used to constrain future parton distribution function fits. 
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PACS numbers: 13.38.Be,13.85.Qk,14.60.Ef,14.70.Fm 



A measurement of the boson rapidity (yw) dis- 
tributions in pp collisions provides valuable information 
about the parton distribution functions (PDFs) of the u 
and d quarks in the proton. W bosons at the Fermilab 
Tevatron pp Collider are primarily produced by quark- 
antiquark annihilation. Contributions from valence- 
valence and valence-sea annihilations provide about 85% 
of the cross section with the rest coming from sea-sea 
annihilations. A W + boson is usually produced by the 
interaction between a u quark from the proton and a d 
quark from the antiproton, while a W~ boson is pre- 
dominantly produced by a d quark from the proton and 
a u quark from the antiproton. Since u valence quarks 
carry on average more of the momentum of the proton 
than d valence quarks 0, 0] , W + bosons tend to move in 
the proton direction and W~ bosons in the antiproton 
direction, giving rise to the W boson production charge 
asymmetry. 

The boson production asymmetry provides com- 
plementary information on the PDFs to that from deep 
inelastic scattering experiments. In particular, it con- 
tributes to determining the slope of the d/u quark ratio 
in the region x < 0.3, where x is the momentum fraction 
carried by a parton in the proton 0, 0] . Knowledge of 
the d/u ratio is needed for predicting the transverse mo- 
mentum (pt) spectrum of leptons from W boson decay, 
an important ingredient in the precision measurement of 
the W boson mass. 

The region of phase space in x that can be probed 
depends on the range of the rapidity of the W boson 
and, at leading order, is given by 



Mi 



w 



(1) 



where a^x(2) is the momentum fraction carried by the u(d) 
quark, +(— )yw is the rapidity of the positive (negative) 
W boson, and y/s is the center of mass energy. For this 
analysis, —2 < yyy < 2 and y^i = 1.96 GeV, allowing us 
to probe 0.005 < x < 0.3. 

The W boson rapidity is not directly measurable since 
the longitudinal momentum of the neutrino from its de- 
cay cannot be determined. However, the rapidity distri- 
bution of the charged lepton from the W boson decay 
reflects this asymmetry; here we use muons for this pur- 
pose. The muon asymmetry is a convolution of the W 
boson production asymmetry and the asymmetry from 
the (V — A) nature of the W boson decay. Since the 
(V — A) interaction is well understood, the muon charge 
asymmetry can be used to probe the PDFs. 

The muon charge asymmetry as a function of the muon 
rapidity is defined as 



A(y) = 



dy 



(2) 



where ^(A i± ) is the differential cross section for the 
boson decay muons as a function of muon rapidity y. In 
this analysis, we measure the muon charge asymmetry as 
a function of pseudorapidity 77 = — ln[tan(f?/2)], where 9 
is the polar angle with respect to the proton beam. 

Allowing for acceptance and efficiency differences be- 
tween positively and negatively charged muons, the muon 
charge asymmetry can be written as 



A(ri) 



(3) 



dy 



where 77 is the pseudorapidity of the muon, N^± (rf) is the 
background-corrected number of muons in pseudorapid- 
ity bin 77, and k(rf) = t + (jf) / e~ (jf) is the acceptance and 
efficiency ratio between the positive and negative muons. 

The lepton charge asymmetry from W boson decay 
was measured by the CDF collaboration in the electron 
and muon channels during Run I of the Tevatron Col- 
lider 0, Q and in the electron channel using Run II data 
0] . The measurement described here is based on a larger 
data sample than these analyses used. In addition, the 
muon channel benefits from a much lower charge misiden- 
tification probability than the electron channel has. 

For this analysis, we used 0.3 fb _1 of data collected 
using the DO detector at a center-of-mass energy of 
1.96 TeV at the Tevatron Collider. The DO detector is 
described in detail in Ref. |8|; here we provide a brief 
description. The detector consists of a central tracking 
system, calorimeters, and a muon detector. The central 
tracking system comprises a silicon microstrip tracker 
(SMT) and a central fiber tracker (CFT), both located 
within a 2 T solenoidal magnet. The SMT has a six- 
barrel longitudinal structure interspersed with sixteen ra- 
dial disks. The CFT has eight coaxial barrels, each sup- 
porting two doublets of overlapping scintillating fibers, 
measuring along the axial direction and at stereo angles 
of ±3°. The liquid- argon/uranium calorimeter is in three 
sections, a central section covering < 1 and two end 
caps extending coverage to \rf\ ~ 4. The calorimeter is 
surrounded by the muon detector which consists of three 
layers of scintillators and drift tubes, one layer in front of 
a 1.8 T iron toroidal magnet and two layers outside the 
magnet. Coverage for muons is partially compromised by 
the calorimeter support structure at the bottom of the 
detector (\rj\ < 1.25 and 4.25 < (j> < 5.15, where <j> is the 
azimuthal angle). 

W — > p,v events were collected using two single-muon 
triggers. An unprescaled trigger, the "wide" trigger, cov- 
ered the region \r]\ < 1.5, and a prescaled trigger, the 
"all" trigger, covered \rj\ < 2. Both triggers required hits 
in the muon system consistent in location and timing 
with a muon originating from the interaction region at 
the first trigger level, a reconstructed track in the muon 
system with px > 3 GeV/c at the second trigger level, 
and a track in the central tracking system with pt > 10 
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GeV/c at the third trigger level. Offline, muon candi- 
dates were required to lie within the geometrical accep- 
tance of the detector, to have reconstructed track seg- 
ments in at least two layers of the muon system, and 
to be matched to a track in the central tracking system 
with transverse momentum pt > 20 GeV/c. To ensure 
well-measured tracks, the matching track was required to 
have at least one hit in the SMT, at least nine hits in the 
CFT, a track fit x 2 /dof < 3.3, and a distance of closest 
approach to the event vertex |dca| < O.Off cm. Muon 
candidates were also required to be isolated in both the 
central tracking system and the calorimeter: the sum of 
the transverse momenta of the tracks in a cone of ra- 
dius 0.5 in r\ — <fi space around the muon track had to 
be less than 2.5 GeV/c, and the total transverse energy 
in a hollow cone of inner radius 0.1 and outer radius 0.4 
around the extrapolated muon position in the calorime- 
ter had to be less than 2.5 GeV. Events with \t)\ < 1.4 
(1.4 < \t)\ < 2) were required to satisfy the wide (all) trig- 
ger. W boson candidates were further selected by requir- 
ing the missing transverse energy ~$ T , determined by the 
vector sum of the transverse components of the energy 
deposited in the calorimeter and the px of the muon, to 
be greater than 20 GeV and the transverse mass Mt > 40 
GeV/c 2 . The analysis was done in muon pseudorapidity 
bins of width 0.2; each bin was treated independently. 
The pseudorapidity resolution of the tracking system is 
approximately 0.01 for \r}\ < 1.7 and is not expected to 
worsen significantly up to \t}\ = 2. 

Events characteristic of the Z — > fifi background were 
removed using two criteria: (i) all events with a second 
muon were rejected, except those within | A0| < 0.1 of the 
selected muon to avoid vetoing events containing multiple 
muons reconstructed from a single real muon at detector 
boundaries, and (ii) all events containing a second track 
with | A(f>\ > 2.1 from the selected muon, pr > 20 GeV/c, 
and |dca| < 0.011 cm were rejected. Cosmic ray muons 
were rejected using muon system scintillator timing cuts 
and by the dca requirement. A total of 189697 W — ► px> 
candidates was selected. 

The asymmetry measurement is sensitive to the 
misidentification of the charge of the muon. A positive 
muon misidentified as a negative muon would not only 
add to the number of negative muons, but would also 
reduce the number of positive muons, and vice versa, 
diluting the true asymmetry. The charge misidentifica- 
tion probability was estimated using a dimuon data sam- 
ple in which events were required to satisfy one of the 
two single-muon triggers described above. Events con- 
taining two muons satisfying all of the above conditions 
(except the second muon veto) and with a dimuon invari- 
ant mass > 40 GeV/c 2 were selected to form a Z 
boson sample. The event sample contained 9958 events 
and only one same-sign dimuon event, giving an average 
charge misidentification probability of (S.O^W) x 10~ 5 
over — 2 < rj < 2. Removing the dimuon invariant mass 
cut or lowering the muon pt cut to 15 GeV/c did not 
change the result. The probability was verified using an 



independently-triggered dimuon sample in which events 
were collected using one of a set of dimuon triggers with 
no track requirements. Out of 19284 dimuon events se- 
lected as described above, two had same-sign muons. We 
also measured the probability using a jsample of W — ► p,v 
events generated with pythia v6.2 [j| and CTEQ6.1M 
PDFs [l0| and passed through the full DO detector simu- 
lation based on GEANT [ll[ ; it is approximately the same 
as that determined using the data. Therefore, charge 
misidentification is not expected to influence the final 
asymmetry measurement. To determine the systematic 
uncertainty on the asymmetry due to charge misiden- 
tification, we used Poisson uncertainties based on the 
number of muon tracks in each r\ bin from the single- 
muon-triggered Z boson sample and varied the asymme- 
try accordingly. The largest uncertainty due to charge 
misidentification is in the range 1.8 < r\ < 2.0, where it 
is 0.001. 

Ideally, the acceptances and efficiencies would be the 
same for all muons as a function of rj and independent 
of charge and pr, leading to k(rf) = 1. To reduce charge 
effects due to detector asymmetries, the directions of the 
magnetic fields in the solenoidal and toroidal magnets 
were regularly reversed. Approximately 51.1% of the se- 
lected W boson sample was collected with the solenoid 
at forward polarity, with 48.9% at reverse polarity. For 
the toroid, 50.7% (49.3%) of the selected sample was col- 
lected with forward (reverse) toroid polarity, respectively. 
In addition, we studied the trigger efficiency for muons at 
the first two trigger levels, the trigger efficiency for tracks 
at the third trigger level, the offline muon reconstruction 
efficiency, the offline tracking efficiency, and the isolation 
efficiency. The four trigger and reconstruction efficien- 
cies are discussed together while the isolation efficiency 
is discussed separately in the following text. 

All efficiencies were measured using the tag and probe 
method on a sample of dimuon events collected using one 
of the single-muon triggers. To select primarily Z boson 
decays, events were required to have M MM > 40 GeV/c 2 . 
First a tag muon was chosen as a track-matched, iso- 
lated muon satisfying all of the selection conditions. The 
probe was another track or muon whose selection criteria 
depended on the efficiency being studied. All efficiencies 
were checked as functions of pr, charge, and r\. No de- 
pendence on pt or charge was observed for the four trig- 
ger and reconstruction efficiencies. The four efficiencies 
were multiplied to determine the combined efficiencies for 
positive and negative muons as functions of rj. Figure [1] 
shows the combined efficiencies by charge and the ratio 
of these efficiencies as functions of rj. The ratio was fit 
to a constant value of 0.99 ± 0.01, with a x 2 /dof of 0.71. 

The isolation efficiency was also measured using the 
tag and probe method on a sample of dimuon events. 
One muon candidate was required to satisfy all muon and 
track selection requirements, while the other had to sat- 
isfy all requirements except that it was not required to be 
isolated in either the tracker or the calorimeter. The frac- 
tion of isolated probe tracks with M MAI > 40 GeV/c 2 gives 
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FIG. 1: (a) The combined trigger and reconstruction effi- 
ciency distributions by charge as a function of r\ and (b) the 
ratio of these efficiencies fit to a constant, 0.99 ± 0.01. 
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FIG. 2: The isolation efficiency as a function of (a) r\ (sep- 
arated by muon charge) and (b) pr (for all muons) fit to a 
constant 0.921 ± 0.002. 



the efficiency. The isolation efficiency is shown in Fig. [2] 
as a function of px and of r\. The efficiency is constant 
for, and consistent between, both charges over the full r\ 
range. The efficiency as a function of px shows a slight 
Pt dependence. We chose, however, to fit the efficiency 
to a constant value of 0.921 ± 0.002 with a x 2 /dof = 5.8. 
To account for the high x 2 /dof, the uncertainty was de- 
termined from the isolation efficiency distribution itself; 
the rms of this distribution is 0.022, and this value was 
used as the isolation efficiency uncertainty. 



TABLE I: Estimated backgrounds in the W boson sample by 
trigger. Uncertainties are statistical only. 



Background 


Wide trigger 


All trigger 


Z —* [ifi 


(4.31 ± 0.05)% 


(4.39 ± 0.11)% 


Z — > TT 


(0.19 ± 0.01)% 


(0.20 ± 0.02)% 


W ->TV 


(2.32 ± 0.02)% 


(2.43 ± 0.08)% 


Multijet events 


(2.77 ± 0.04)% 


(2.76 ± 0.09)% 



The largest source of contamination in the data sample 
comes from electroweak processes, Z — > fifi and W — > tv 
and Z — ► tt where r — > [ivv. Muons from these elec- 
troweak decays exhibit charge asymmetries which di- 
lute the true asymmetry. These asymmetries are ac- 
counted for by subtracting the background bin-by-bin 
in r\. The electroweak background was estimated us- 
ing Monte Carlo samples generated with PYTHIA v6.2 
and CTEQ6.1M PDFs and a parameterized description 
of the DO detector !12|. For each of the three processes, 
separate samples were generated for the two triggers. 

The background due to semi-leptonic decays and 
punch-through in multijet events was estimated using 
the data. The isolation criteria remove events containing 
muons within jets, and they were used as the discrimi- 
nator to determine this background. Using a sample of 
events passing all selection criteria except that on the 
transverse mass and requiring Jfi T < 10 GeV, we mea- 
sured the probability for multijet events to satisfy the 
isolation criteria. This probability shows no dependence 
on muon pseudorapidity. 

Table U shows the overall contribution of each of the 
four backgrounds. To determine the number of events for 
each background, the contributions from the three elec- 
troweak processes were added to the number of events 
expected from W — > [iv Monte Carlo events produced us- 
ing PYTHIA and the parameterized detector description, 
and this sum was normalized to the number of events 
in the data less the estimated multijet background. The 
overall normalization was done for |?y| < 1.5 for the wide 
trigger and \r]\ < 2 for the all trigger, while for the final 
result, it was done independently in each r\ bin. This is 
the only use of signal Monte Carlo events in the analysis. 

The muon charge asymmetry was determined sepa- 
rately for each bin in rj and is shown in Fig. [3] Also 
shown are the asymmetry determined using the resbos 
[lH event generator with QCD resummation and PHO- 
TOS NLO QED corrections in the final state with the 
CTEQ6.1M PDFs^with the forty CTEQ6.1M PDF un- 
certainty sets [13, El, and with the MRST04 NLO PDFs 
[lt| . The next-to-next-to- leading-order calculation of the 
W boson rapidity distribution [l7[ at the Tevatron is 
very similar to the next-to- leading-order calculation. The 
W boson asymmetry distribution is not very sensitive 
to QCD corrections, and calculations at leading-, next- 
to-leading-, and next-to- next-to-leading-order are nearly 
indistinguishable [l7| . 

Systematic uncertainties taken into account are those 
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TABLE II: The measured muon asymmetry in bins of pseu- 
dorapidity, calculated using Eq. [3] The asymmetry values 
are the averages within each pseudorapidity bin. The first 
uncertainty is statistical; the second is systematic. 

Pseudorapidity range Muon asymmetry 
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n a 
— U.4 — 


— U.z 


— U.U4y ztz U.U11 ztz U.UUO 


-0.2- 


0.0 


-0.011 zb 0.010 =b 0.005 


0.0- 


0.2 


0.028 ± 0.011 ± 0.005 


0.2- 


0.4 


0.050 =b 0.011 =b 0.005 


0.4- 


0.6 


0.071 =b 0.011 =b 0.005 


0.6- 


0.8 


0.120 zb 0.011 zb 0.005 


0.8- 


1.0 


0.122 zb 0.012 zb 0.005 


1.0- 


1.2 


0.127 =b 0.011 =b 0.005 


1.2- 


1.4 


0.107 =b 0.009 =b 0.005 


1.4- 


1.6 


0.065 zb 0.025 =b 0.007 


1.6- 


1.8 


0.042 zb 0.036 =b 0.005 


1.8- 


2.0 


-0.102 zb 0.087 zb 0.005 
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FIG. 3: The muon charge asymmetry distribution. The hor- 
izontal bars show the statistical uncertainty and the full ver- 
tical lines show the total uncertainty on each point. The 
shaded (yellow) band is the envelope determined using the 
forty CTEQ6.1M PDF uncertainty sets, the solid (red) line 
is the CTEQ6.1M central value, and the dotted (blue) line is 
the charge asymmetry determined using the MRST04 NLO 
PDFs. All three were determined using RESBOS and PHOTOS 
(color online). 



on the ratio of the efficiencies for positively and nega- 
tively charged muons [^(ry)], the efficiency of the isola- 
tion criteria, the charge misidentification, the probabil- 
ity for multijet events to satisfy the isolation criteria, the 
hadronic energy scale for the detector (needed to calcu- 
late $ T ), and the parameterization of the muon energy 
loss in the calorimeter. Each contribution was varied by 
±lcr, and the asymmetry was recalculated. The changes 
in the measured values were added in quadrature to de- 
termine the overall systematic uncertainty. The mea- 
sured asymmetry and the statistical and systematic un- 
certainties are listed in Table HT1 for each r\ bin. The ef- 
ficiency ratio uncertainty, which is based on the number 
of dimuon events in the single- muon-triggered samples, 
dominates and is approximately equal to 0.005 in all r\ 
bins. 

By CP invariance, the asymmetries at ±77 have oppo- 
site signs and equal magnitudes, allowing the asymmetry 
distribution to be "folded" to decrease the statistical un- 
certainty. The folded asymmetry distribution was found 
by combining the numbers of events in each ±77 bin and 
redetermining the background and systematic uncertain- 
ties as described above. The folded distribution is shown 
in Fig. 3] with the measured values of the asymmetry and 
uncertainties given in Table [TTTT For 0.7 < |?7| < 1.3, our 
experimental uncertainties are smaller than the uncer- 
tainty given by the CTEQ uncertainty sets. Only at the 
extremes of our measurement, in the 0.0-0.2 and 1.8-2.0 
muon pseudorapidity bins, are our uncertainties larger 
than the CTEQ uncertainty. Between these region, the 
uncertainties are comparable. 

We have measured the charge asymmetry of muons 




Muon pseudorapidity 

FIG. 4: The folded muon charge asymmetry distribution. The 
horizonal bars show the statistical uncertainty and the full 
vertical lines show the total uncertainty on each point. The 
shaded (yellow) band is the envelope determined using the 
forty CTEQ6.1M PDF uncertainty sets, the solid (red) line 
is the CTEQ6.1M central value, and the dotted (blue) line is 
the charge asymmetry determined using the MRST04 NLO 
PDFs. All three were determined using RESBOS and PHOTOS 
(color online). 



from W boson decay using 0.3 fb _1 of data. Our results 
can already improve constraints on the PDFs. In the fu- 
ture, measurement of the W boson asymmetry will have 
a significant impact on PDF determination as present 
uncertainties are dominated by statistics. 

We thank the staffs at Fermilab and collaborating in- 
stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
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TABLE III: The folded muon asymmetry in bins of pseudo- 
rapidity. The asymmetry values are the averages within each 
pseudorapidity bin. The first uncertainty is statistical; the 
second is systematic. 

Pseudorapidity range Muon asymmetry 



0.0- 


0.2 


0.019 ±0.008 ±0.005 


0.2- 


0.4 


0.050 ±0.008 ±0.005 


0.4- 


0.6 


0.081 ±0.008 ±0.005 


0.6- 


0.8 


0.126 ±0.008 ±0.005 


0.8- 


1.0 


0.121 ±0.008 ±0.005 


1.0- 


1.2 


0.133 ±0.008 ±0.005 


1.2- 


1.4 


0.124 ±0.006 ±0.005 


1.4- 


1.6 


0.085 ±0.017 ±0.006 


1.6- 


1.8 


0.031 ±0.026 ±0.005 


1.8- 


2.0 


-0.006 ±0.062 ±0.005 
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